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This paper presents a mathematical analysis of two-dimensional steady-state heat 
conduction in a solid of rectangular cross section, two parallel surfaces of which are exposed 
to separate temperature environments with respective surface heat transfer coefficients. 
For the other two parallel surfaces, the temperature environment and surface heat transfer 
coefficients are assumed to vary as a function of position on these surfaces. A particular 
application of this analysis has been to determine the temperature distribution in a trans- 
verse cross section of exposed concrete columns. 



1. Introduction 

Modern architectural treatment of exterior columns of high-rise buildings involves the 
use of columns exposed to temperature extremes on the outside face and relatively constant 
temperatures on the interior. The temperature extremes give rise to dimensional changes in 
the column which can best be determined from temperature distribution in the column. In 
connection with a request from the Federal Housing Administration, an investigation was 
made for determining the temperature distributions in an exposed exterior column. 

This paper presents a mathematical analysis of two-dimensional steady-state heat con- 
duction in an exposed exterior column which may be represented by a solid of rectangular 
cross section, two parallel surfaces of which are exposed to separate temperature environments 
with respective surface heat transfer coefficients. For the other two parallel surfaces, the 
temperature environment and surface heat transfer coefficient are assumed to vary as a function 
of transverse position on these surfaces. 

An analysis for general application is presented, along with the statement and numerical 
results for four specific problems. 

2. Analysis for General Problem 

Because the heat-transfer problem in a rectangular solid may have a more genera] applica- 
tion than that cited in the previous section, the following analysis is based on boundary con- 
ditions, arbitrary in nature, which may be adapted to specific applications. Section 3 presents 
the assumed boundary conditions for heat transfer in an exposed exterior column. 

For steady-state, two-dimensional heat conduction in a rectangular homogeneous solid 
with no flow of heat in the direction of the third dimension, it is assumed that the rate of heat 
flow at each of the four boundary surfaces is proportional to the temperature difference be- 
tween the temperature of the surface and that of the ambient adjacent to the surface. The 
proportionality factor may be prescribed as a function of the distance along the surface. A 
•convenient placement of coordinate axes for the type of problem to be considered is as shown 
in figure 1, where the boundary conditions at the four surfaces, x= ±Z, and y= ±a, are shown 
adjacent to the surfaces. 
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Figure 1. Transverse section of rectangular solid, 
showing boundary conditions on x=±l and 
y=±a. 



The temperature potentials V and V(x) are the temperature differences between the 
temperatures of the ambient adjacent to various surface positions and that of the ambient 
adjacent to x=l; the temperature potential v is the temperature difference between a tempera- 
ture in the solid and that of the ambient adjacent to x=l. The proportionality constants 
H a and H e represent the coefficients of heat transfer between the surface and the ambient at 
the surfaces x=—l, and x=l, respectively. For the purposes of this paper the variable heat 
transfer coefficient H(x) and temperature potential V(x) are considered to be piece-wise con- 
tinuous on —l<x<l, where H{x)=H t {x) and V(x) = Vi(x) are continuous on Xi_i<x<Xi and 
i=l, 2, 3, . . . m; # =— Z, and x m —l. 

The partial differential equation for steady-state heat conduction in a homogeneous solid 
in a rectangular coordinate system is 



b 2 v . d 2 fl__ n 
dx 2 dy 2 

A solution of (1) satisfying the boundary conditions at x=l and x=—l (fig. 1) is 

v _hi+hih 2 (l—x/l) , _ A n R(p n x/l) cosh p n yfl 



V h 1 +h 2 +2h 1 h 2 



%=i 



cosh (3 n a/l 



(i) 



(2) 



where due to the similar conditions on y= ±a, it is assumed that i-=0on the center line y=0. 
Symbols used in (2) are: 

K= thermal conductivity of the solid 

h 2 =H e l/K 
A n = coefficients to be determined from boundary conditions on y=a 

R(p n x/l) = p n cos p n (l+xll) + h l sin p n (l+x/l) 

where /3 n are positive roots of (& 2 —liih 2 ) sin 2/3— j8(Ai+A 2 ) cos 2/3=0. 
Substitution of boundary conditions on y=a gives m equations 

Multiplying each equation by R(p k x/l) and integrating from x= — l to x=Z, yields 
S I iMU tanh ^a/ZJ_^(^x/Z)i?(^/Z)cZx+-^g ^J H(x)R((3 n x/l)R(p k x/l)dx I 
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Using the orthogonal and other related properties of R(p n x/l) as derived in the appendix, (4) 
may be put in the form of simultaneous equations 



where 



£>nA n -j-2__, C n>k A k — D n , n — 1, 2, . . . 

k=l 
kj*n 

C l 1 m ( % x i 

B n =p n tanh n a/l\ R 2 (i3 n x/l)dx+~J^ H(x)R 2 (p n x/l)dx 

0,,*=4f) P H(x)R(p n xll)R(i3 k x/l)dx 



B -ls£,[ 



V(x) h 1 +hMl~x/l) 



f]x« 



x)R(p n x/l)dx. 



V h 1 +h 2 +2h l i 
Of particular interest is the mean temperature in the solid, # OT , given by 

v m h x (l+h 2 ) . I ^A n tB,i)hp n a/l[p n sm2p n +h 1 (l--cos2p n )] 



V hi-\-h 2 J r2h l h 2 2a » 
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3. Temperature Distribution in Exposed Concrete Columns 



(5) 



(6) 



An interesting specific case for solution of (2) is to consider a transverse section of a rec- 
tangular concrete column. Here a portion of the boundary surface is exposed to outside weather 
conditions or design temperature, another portion of the surface is exposed to the ambient 
temperature maintained within the building, and another portion of the surface is in contact 
with the walls of the building. A cross section of the column, and an abutting wall in the region 
b<^x<^d, is shown in the lower portion of figure 2. The x=l face of the column is assumed to be 



Figure 2. Transverse section of rectangular concrete 
column with abutting wall separating the inside and 
outside environment. 

Outside face is located at x =1. 
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the outside face exposed to the outside design temperature which is assumed to be the zero or 
datum plane temperature from which all other temperature potentials are reckoned; the face 
x=—l is assumed to be exposed to the indoor temperature potential V. 

The upper portion of figure 2 shows the assumed variation of H(x) and V(x) along the 
surfaces y=±a. H(x) is assumed constant at H a in the region — l<x<b; decreases linearly 
from H a to zero in going from b to c; increases linearly from zero to H e in going from c to d; 
and is constant at H e in the region d<x<l. V(x) is assumed constant at V from —I to c; 
and is zero from c to I. In the region of the abutting wall, 6<x<rf, the heat transfer between 
the wall and the column was assumed to be dependent on the inside and outside temperatures 
in their region of influence and the surface coefficients decreasing linearly to zero at x=c. 



dy 



(SJllf ) (*-^) b<x<c, y=a 
~ K %= H *(jES) v c<x<d iV =a. 



Expressions for D ni B n , and C n< k necessary for substitution in (5) are given in the appendix 
(A13, Al4, A15). The solution of (5) for A n and its substitution in (2) gives the temperature 
distribution for values of x and y. 

Numerical solutions were obtained, using a digital computer, for a 36 in.X14 in. concrete 
column for four positions of the abutting wall, based on a 100 deg F temperature difference 
between the inside and outside. Numerical values used, and assumed as constants, were 

Thermal conductivity, If =1.0 Btu/hr ft deg F 
Inside heat transfer coefficient, H a = 0.5 Btu/hr ft 2 deg F 
Outside heat transfer coefficient, H e =6.0 Btu/hr ft 2 deg F 
Width of column, 2a=14 in. = 7/6 ft 
Thickness of column, 2/=36 in. =3 ft 
Thickness of wall, d—b=8 in. =2/3 ft 
Also it was assumed that d—c=c—b. 

The inside surfaces of the column were assumed to be separated from the inside of the 
building by suitable interior finish, giving a nominal value for H a of 0.5 Btu/hr ft 2 deg F. The 
value i? e =6.0 Btu/hr ft 2 deg F was based on an outside wind velocity of 15 mph. 
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Figure 3. Temperature potential distributions in 36-in. by 11^-in. concrete columns for four positions of the 
abutting walls, based on a 100-deg temperature difference between the inside and outside and the constants 
cited in the text. 
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Figure 3 shows temperature potential isotherms plotted within the column cross section 
for four positions of the abutting wall. Mean temperatures were computed from (6). Fig- 
ure 4 shows mean temperatures over the width of the columns plotted against distance meas- 
ured from the inside surface of the column. Positions of the abutting wall are shown for 
clarity. 

Temperatures in degrees Fahrenheit can be expressed in terms of the temperature potentials, 
v, of figure 3 by means of the relation 

t, °F = t e +(v/lOO)(t a -t e ) 

where t a is the inside air temperature, and t e is the outside air temperature, in degrees 
Fahrenheit. 

4. Discussion 

The analysis of section 2 assumes similar boundary conditions on the surfaces, y=±a. 
For the consideration of a more general case where these conditions are not similar, it is con- 
venient to make a translation of the x-axis of figure 1 to the lower side and rewrite the sum- 
mation term in (2) 

^A n COSh PnVll+En COsh $ n (^J^)] R(Ml) 

cosh (3 n 2a/l 

Substitution of boundary conditions at y=0 and y=2a will give two sets of simultaneous 
equations similar to (5), both containing the coefficients A n and E n . Solution for A n and 
E n may be obtained by iteration; for example, let E n =0 in one set and solve for A n , substi- 
tute A n in the second set and solve for E n , which is substituted in the first set. This process 
may be continued until A n and E n are invariant. 



Figure 4. Mean temperatures over the width of the 
columns against distance measured from the inside 
surface of the columns. 

Positions of the abutting walls are shown. 
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The assumption was made that a sufficiently good approximation to the solution of the 
infinite set of equations (5) in the infinitely many unknowns A n could be obtained by solving 
the first forty equations for A n , n=l, 2, . . ., 40. Solution of the first twenty equations for 
A n , 7i=l, 2. . ., 20 gave results which differed from that for forty A n by less than 0.1 °F. It 
therefore appeared that convergence of the approximating sequence of A n 's to the true solution 
of (5) was sufficiently rapid to justify our assumption. 

A computer program in Fortran IV language is available for the problem as set forth in 
section 3. Data input to this program allows any possible variation in the numerical values 
for the thermal conductivity, surface coefficients, and dimensions. 

4. Appendix. Orthogonal and Other Related Properties of R(fi n x/l) 
To satisfy the boundary condition in (2) at x=l, fi n must be the roots of 

(P 2 —Kih) sin 2/3-/3(^+^2) cos 2/3=0. (Al) 

Without proof in this paper, the roots fi n will be assumed real and positive. Proof of the real 
nature of these roots is similar to that for roots of an expression given in Carslaw and Jaeger. 1 
It then becomes expedient to develop the orthogonal properties of 

R(i3 n x/l) = (3 n cos p n (l+x/l)+h 1 sin p n (l+x/l) (A2) 

where — Z<!x<CZ. Several expressions useful for later developments are cited below: 

sin 2 og Mhi + h,)' (A) 

sin 2P- { p +mj32+m {A3) 



(/3 2 -W 2 



cos 2/3— ! 2 . , 2 \ i «2 j m (A4) 



E\fi x /l)=^~+^—^ cos 2/3(1 +x/l)+h£ sin 2/3(1 +z/Z). (A5) 

R(Ml)R(MD= Ll{x) ^ L2{x) (A6) 

L x (x) = {PA—hX) cos«(l+a;/Z)+M sin <t>(l+x/l) 
L 2 (x) = (PA+h 2 i)cosf(l+x!l)-h 1 +sin+(l+x/l) 

M a (x)= (R(p n x/l)dx=j- [ft, sin /3,(l+x/Z)-A! cos ft^l+z/Z)] (A7) 

No(x)=jB\0 n x/l)dx= x{ ^ h2l) + lifi ljf l) sin 2ft i (l+x/Z)-^ cos 2/3,(1+*/*) (A8) 

P (x)=JL l (x)dx= l - [(ftA-A?) sin ^(l+x/D-h^ cos <t>(l+x/l)] (A9) 

Qo(x)= JL 2 (x)dx=j [(ftA+AJ) emtf.l+x/l)+W cos*(l+x/Z)]. (A10) 



Carslaw, H. S., and Jaeger, J. C, Conduction of Heat in Solids, 2d ed., p. 114 (Oxford University Press, London, 1959). 
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By use of (A3) and (A4), it can be shown that P (l) + Q (l)=Po(—l) + Qo(—l)=0] i.e., the 
integral R((l n xll)R(p k x/l)dx=0, showing that R(fi n x/l) and R(p k x/l) are orthogonal to each 
other in the interval —l<Cx<^l. Also, using (A3) and (A4), 

No{l) = l ($l+h\){Pl+ht+h 2 ) 



2 (fil+hl) 

Assuming H t (x) and Vi(x) to be polynomials in x, then the recurrence 

M p (x)= \x p R(px/l)dx=xM p . 1 (x)- \M P ^(x)dx (All) 

and similar relationships for N p (x), P p {x) and Q p (x) may be established for determining B n > 
D ny and C n>k in (5). Assuming a resulting polynomial of the form S-\-Tx-\- Ux 2 in the solution 
for D ny the indefinite integral becomes 

((S+Tx+Ux 2 )R(Ml)dx=(s+Tx+Ux 2 -^) (A12) 

./ \ Pn J Pn 

The solution for T) n for the specific case illustrated in figure 2 takes the form 



D n =- t 



h x h 2 



T S ( X {S i + T i x+U i x 2 )R{fi n xll)dx 

i/h i=l JXi-i 



h l +h 2 +2K 
where 

S 1 =l+h 1 2\=M Z7i=0 

S 2 =cS 1 l(c-b) T 2 ^{c r I\-S x )l{c-b) U 2 =-T 1 /(c-b) 

S*=-cSJ(d-c) T 3 =(S,-cT 4 )/(d-c) U 3 =TJ(d-c) 

St=-(l+h 2 ) T A =h 2 /l U 4 =0. 

The function in the region — /<C:rO is a linear, piecewise continuous function for which 
the integrals go to zero at the integration limits x=—l and x=l, and elsewhere at the other 
limits only the coefficients of l 2 /fil in (A12) need be considered, 

+(T l +T i -2U i b)R(f} n b/l)-(T i -T 3 -2U 3 d)R(J3 n d/l) 

+ {T r -T a +2c({7 t -U t ) } HiMDl (A13) 
By a similar analysis 

g» = (ft. tanh p n a/l+h2Wn+h 2 l )(Pl+hi+h 2 ) + (0 n tanh g.o^+fti)(|8.+^)03.+ft?+Ai) 
* " 2(/3*+A|) 

+ iS [(&) { R2 ^dJl)-R 2 (Ml) } -^y {R*(Ml)-BW/l) }] 



p t (c+6) h(d+e) l (AH) 



4 

^ ^iW-^M , L 2 {d)-L 2 {c) -\ , A,? r ^i(«)-A(i) , X a (c)-Z 2 (6) "l , . . .. 
~2(d-c)L 4> 2 * 2 _T2(c-6)L <t> 2 ~ + ~ ^ J' (A 5j 

(Paper 69C2-194) 



151 



